Non-Gaussianity in WMAP Data Due to the Correlation of CMB Lensing Potential 

with Secondary Anisotropies 
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We measure the skewness power spectrum of the Cosmic Microwave Background (CMB) 
anisotropies optimized for a detection of the secondary bispectrum generated by the correlation 
of the CMB lensing potential with integrated Sachs- Wolfe effect and the Sunyaev-Zel'dovich ef- 
fect. The covariance of our measurements is generated by Monte-Carlo simulations of Gaussian 
CMB fields with noise properties consistent with Wilkinson Microwave Anisotropy Probe (WMAP) 
5-year data. When interpreting multi-frequency measurements we also take into account the con- 
fusion resulting from unresolved radio point sources. We analyze Q, V and W-band WMAP 5-year 
raw and foreground-cleaned maps using the KQ75 mask out to Zmax = 600. We find no significant 
evidence for a non-zero non-Gaussian signal from the lensing-secondary correlation in all three bands 
and we constrain the overall amplitude of the cross power spectrum between CMB lensing potential 
and the sum of SZ and ISW fluctuations to be 0.42 ± 0.86 and 1.19 ± 0.86 in combined V and 
W-band raw and foreground-cleaned maps provided by the WMAP team, respectively. The point 
source amplitude at the bispectrum level measured with this skewness power spectrum is higher 
than previous measurements of point source non-Gaussianity. We also consider an analysis where 
we also account for the primordial non-Gaussianity in addition to lensing-secondary bispectrum and 
point sources. The focus of this paper is on secondary anisotropies. Consequently the estimator is 
not optimised for primordial non-Gaussianity and the limit we find on local non-Gaussianity from 
the foreground-cleaned V-l-W maps is /nl = — 13 ± 62, when marginalized over point sources and 
lensing-ISW/SZ constributions to the total bispectrum. 
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PACS numbers: 98.80.-k 95.85. Sz, 98.70.Vc, 98.80.Cq 

I. INTRODUCTION 

The all-sky, multi-frequency WMAP maps of Cosmic 
Microwave Background (CMB) anisotropics have pre- 
sented cosmologists with an opportunity to test the cos- 
mological scenario of structure formation at an unprece- 
dented accuracy. The results on the CMB temperature 
and polarization angular power spectra are in very good 
agreement with the expectations of the standard cosmo- 
logical model of structure formation based on primor- 
dial, adiabatic and scale invariant perturbations, and 
cold dark matter ||]. 

In addition to measuring the angular power spectrum 
and cosmological parameter estimates from it |3| , WMAP 
CMB maps are now routinely used to constrain sta- 
tistical properties of the CMB beyond the simple two- 
point angular correlation function. These studies in- 
clude tests of cosmological isotrop y JM p i H i, topol- 
ogy 01 ^ > ^-nd non-Gaussianity [|| ]l0[~pTr [l^ , among 
other tests. In the standard cosmological model, pri- 
mordial CMB anisotropies are supposed to be Gaus- 
sian, however non-Gaussianities may be present in the 
observed CMB maps through a combination of primor- 
dial non-Gaussianity of density perturbations generated 
during inflation ^ [l^, the imprint of non- 



linear growth of structures as probed by secondary effects 
p8| , |l9|, and mode-coupling effects by secondary sources 
of temperature fluctuations such as gravitational lensing 
of the CMB 1^, |2l|. 

The detection of these non-Gaussian features will not 
only provide additional useful information on the pa- 
rameters of the standard cosmological model but also 
allow independent tests on constraining the amplitude 
of primordial non-Gaussianity due to non-standard ini- 
tial conditions and, ultimately, inflation after account- 
ing for secondary non-Gaussian signals generated since 
last scattering. Several recent studies have made use of 
the bispectrum for the study of non-Gaussianities ^ . 
This quantity involves a three-point correlation function 
in Fourier or multipole space. The configuration depen- 
dence of the bispectrum B{ki, k2, k^) with lengths (fci, 
^2, fca) that form a triangle in Fourier space can be used 
to separate various mechanisms for non-Gaussianities, 
depending on the effectiveness of the estimator used. 

In most CMB non-Gaussianity studies |l^, ^ the 
estimator employed involves a measurement that com- 
presses all information of the bispectrum to a single 
number called the cross-skewness computed with two 
weighted maps. Such a drastic compression limits the 
ability to study the angular dependence of the non- 
Gaussian signal and to separate any confusing fore- 
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grounds from the primordial non-Gaussianity. More re- 
cently, some of us have introduced a new estimator that 
pre serves some angular dependence of the bispectrum 
|2^ . This recently led to a new measurement of the 
primordial non-Gaussianity parameter |^ . 

The skewncss power spectrum is indeed a weighted 
statistic that can be tuned to study a particular form 
of non-Gaussianity, such as what may arise either in the 
early Universe during inflation or late-times during struc- 
ture formation, while retaining information on the nature 
of non-Gaussianity more than the skcwness alone. When 
applied to the CMB data, this allows a way to explore all 
non-Gaussian signals, including those generated by con- 
taminants such as Galactic foregrounds and unresolved 
point sources. 

In this paper we analyze the recent WMAP data 
for the skewness power spectrum associated with cross- 
correlation between lensing and ISW and SZ effects, re- 
spectively. The presence of a measurable signal in this 
secondary non-Gaussianity, especially with the cross- 
correlation of lensing with the SZ effect, was identified 
in 2003 by Goldberg & Spergel @ H. We provide first 
constraints on this signal from WMAP data using Q, V 
and W-band maps both in the "Raw" and foreground 
"Clean" form as provided by the WMAP team publicly. 

After accounting for the confusion from point sources 
generated by the overlap of the point source shot-noise 
bispectrum and the lensing-secondary anisotropy cross- 
correlation bispectrum, we find no significant detection 
of the lensing effect in existing WMAP data. We con- 
strain the overall normalization of the lensing-SZ and 
lensing-ISW angular cross-power spectra to be 0.42 ± 
0.86 and 1.19 ± 0.86 in combined V and W-band raw 
and foreground-cleaned maps provided by the WMAP 
team, respectively. The point source amplitude we de- 
termine from the raw map of Q-band with bgic ~ (67.8 ± 
5.4) X 10~^^ sr^ is higher than the estimate by the WMAP 
team with (6.0 ±1.3) x 10"^ ^K^-sP [§ (the value we 
determine is (13.7± 1.1) x 10"^ fiK^ -sr^ in the same units 
used by the WMAP team) . We find similarly a factor of 
2 increase in the results from the V-band map. In the 
case of clean maps, we find bgyc = (6.2 ± 5.4) x 10^^^ sr^, 
which is smaller than the WMAP team's estimate with 
clean maps for the Q band with (4.3±1.3) x 10~^ ^K^-sr^ 
[H (the value we determine is (1.4 ± 1.1) x 10"'^ ^K^-sr^ 
in the same units used by the WMAP team). We find 
similar differences in the V and W bands as well. It is 
unclear exactly where these differences come from. We 
do not employ the same E-statistic that is optimized for 
point sources as the WMAP team in the present study. 

We also considered the extent to which primordial non- 
Gaussianity confuse the detection of lensing-secondary 
correlations and found that when including /nl hi model 
fits, in addition to point sources, leads to a factor of 
2 degradation in the error of the amplitude of lensing- 
secondary correlation power spectrum. 

This paper is organized as follows: in the next Section, 
we review the measurement theory. We refer the reader 



to Munshi et al. for more details. In Section III 

we present our results and discuss the evidence for the 
secondary non-Gaussianity in WMAP data. 



II. SKEWNESS POWER SPECTRUM 
ESTIMATOR 

If we consider three statistically isotropic fluctua- 
tion fields, say temperature anisotropics but weighted 
with different window functions differently, X(fl), 
Y{il) and Z{(l) described by the multipole moments 
a^^^ , a^^^, fl;^^^, all the information available in the 
three-point correlation function is contained in the angu- 
lar bispectrum B^^^^ defined by a triangle in multipole 
space with lengths of sides (/i, I2, 13) : 
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Since measuring the full bispectrum is challenging, 
many previous measurements have focused mostly on the 
skewness which is collapse of information in the bispec- 
trum in some way to a single number. As discussed in 
Munshi et al. j2^, it is useful to pursue instead the skew- 
ness power spectrum which can be considered as the an- 
gular power spectrum of the correlation of the product 
map X{d)Y{d) and the Z{Cl). In the absence of sky-cut 
and instrumental noise, we can write the skewness power 
spectrum as : 



/ XY z 
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where is the spherical harmonic transform coeffi- 

cient of the field XY . 

It is possible to show that this quantity, in the homo- 
geneity and isotropy assumption, is directly connected 
with the mixed bispectrum associated with these three 
fields according to the relation p3 : 
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where wi-^^i^ is a filter function that needs to be introduced 
in a more general approach and represents the spherical 
transform of the mask. This power spectrum contains 
information about all possible triangular configuration 
when one of the side is fixed at length I. 

We can now relate the X, Y and Z fields introduced 
above to quantities that we are interested in studying. 
We therefore expand the observed CMB temperature 
anisotropics 6T{Q) in terms of the primary anisotropy 
STp, due to lensing of primary 5Tl, and the other secon- 
daries generated by the low-redshift large-scale structure 
STs : 



STin) = 6Tp{Cl) + STl{Ci) + STsin) 



(4) 



3 



Expanding these fields in the Fourier space according 



to 



Im 

jTi(ii) = ^ve(fi)-vT5(fi), 

aSlTnYlmi^) 



(5) 



we have an expression for the cross-correlation power- 
spectra which denotes the coupling of lensing with a spe- 
cific form of secondary CMB anisotropics. Their bispec- 
trum is given by : 
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where {dT)i„i is the anisotropy map expansion to multi- 
pole harmonics ||2^, |2^, With explicit calculations, 
the bispectrum becomes : 
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where Xi^ is the lensing potential and secondary 
anisotropics cross-correlation power spectrum and C/^ is 
the unlensed power spectrum of CMB anisotropics. 



A. Optimised sltew spectrum 

Following the discussion in Munshi et al. we de- 
fine a set of 9 different fields of weighed temperature : 
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where bi is the beam transfer function; Ni is the noise 
power spectrum as obtained from averaging noise maps 
simulations; C; and Ci are the unlensed and the Icnscd 
CMB power spectrum, respectively. 

From these harmonic coefficients, we also construct 9 
sky maps : 
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where i = 1,2, 3. 

The skewncss power spectrum that is weighted to mea- 
sure non-Gaussianity associated with lensing-secondary 
correlation can be written as : 
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The above form is exact for all-sky measurements. To 
account for partial sky coverage due to the Galactic and 
foreground mask and inhomogeneous noise, we also cal- 



culate the linear-order correction terms from Ref. 
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The term above without an averaging is the direct esti- 
mate from data while the averaged corrective terms such 
as cf^^'""^ are obtained by cross-correlating the prod- 
uct of the observed A map with the simulated B and C 
maps and then taking an ensemble average over many 
realizations. The reduction in the sky are due to mask is 
corrected dividing by the observed sky fraction f sky- 
As discussed in Ref. |2^, it is possible to show that 
this quantity is directly related to the bispectrum : 



Cf'' = 
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where Bu-^i^ is the total bispectrum in the data and 
(-^i^i^f ) is the shape of the bispectrum that we have 
employed by weighting the A, B and C maps. This bis- 
pectrum is equal to the form written in equation (fzh , with 
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permutations only restricted to h I2 while I3 is kept 
fixed to Xi.^. 

We assume the total bispectrum present in the data is 
a combination of the both the lensing-secondary bispec- 
trum and contaminations such as point sources. When 
fitting to measurements, we will construct the map C" in 
above by appropriately weighting it with Xi to study the 
cross-correlation of lensing potential with SZ and ISW 
separately. We assume that the bispectrum in the data 
is composed by these two effects with two unknown am- 
plitudes relative to the prediction under the fiducial cos- 
mological model. The comparison between the data and 
the modeled expectation will be used to determine the 
two relative 



III. DATA ANALYSIS 



We summarize our analysis in the following steps: 



A. Data and Simulations 

We have considered the WMAP 5-year Stokes-I raw 
and clean sky maps for the Q, W and V frequency bands, 
obtained from the public lambda website^. We use the 
Healpix's synfast code ||2^ to generate 250 CMB tem- 
perature anisotropy Gaussian maps giving in input the 
WMAP 5-year best-fit CMB anisotropy power spectrum 
with cosmological parameters: Hq = 71.9 km/s/Mpc, 
flbh^ = 0.02273, flch^ = 0.1099, = 0.963 and 
T — 0.087. We require Nside = 512 and a maximum 
multipole equal to 600. 

In the same way, we generate 250 noise maps with noise 
properties consistent with WMAP Q, W and V frequency 
bands : 



O"0 



(13) 



where N{n) is the final noise map obtained from a 
white noise map n(Cl) combined with the WMAP rms 
noise per observation, ao, and the number of obser- 
vations per pixel, A^obs- We extract Nobs from the 
WMAP 5-year Stokes-I sky maps fits files and take 
(To = 2.197,3.133,6.538 mK for the Q, V and W 5-year 
data, respectively, as published on the lambda website 
by the WMAP team. 

We analyze both raw maps as well as foreground- 
cleaned maps provided by the WMAP team. We show 
and tabulate results separately for these two options. 

We use the Healpix anafast code ||2^ and the KQ75 
mask to extract the multipole coefficients for each fre- 
quency band out to Imax = 600 for WMAP maps. 



af^y,, simulated Gaussian maps, a^, and simulated noise 
maps, a/^. The noise spectrum needed for computing 
the denominators in (H) is obtained averaging the simu- 
lated noise spectra over the solid angle and considering 
the sky-cut according to the relation : 



(Ph al M (n) 

AnfskyNobsin) 



(14) 



^ http://lambda.gsfc/nasa.gov 



where fi = Air/Npixei is the solid angle per pixel, M{n) 
is the KQ75 mask and fsky = 0.718 is the corresponding 
observed sky fraction. 

We calculate the lensed and unlensed CMB power spec- 
trum with the public CAMB code |^ using again the 
cosmological parameters from the WMAP 5-year best fit 
model. 

We put everything together to obtain all coefficients 
in (^) and the relative sky maps considering aim = afm 
for data instead, in the case of simulations, we need 
to consider noise and beam contribution to multipoles: 
dim = o,friM + o,iin^ gaussian multipoles are con- 

volved with the frequency-dependent beam transfer func- 
tion hi and added to the noise multipoles. 



B. Skewness power spectrum 

We estimate the C^'"^ according to equation (|ll|) for 
each frequency band and for different lensing-secondary 
anisotropy cross-correlation power spectrum; in particu- 
lar, Xi ISW is the spectrum of cross-correlating lensing 
with Integrate Sachs- Wolfe effect and, in the same 
way, Xi SZ for Sunyaev-Zel'dovich effect jH] (see Fig- 
ure |l|). We calculate these in the fiducial cosmological 
model consistent with WMAP 5-year data and making 
use of the halo model approach to describe the SZ sig- 
nal |l], m. The ISW effect is described through 
the standard linear power spectrum of the potential field 
and the CMB lensing potential is also modeled using the 
linear fluctuations [33], 

In Figures H and |we show all terms of equation (pi]). 
It's evident that linear terms are not significant compared 
to the others coming from data only. However, we are 
using all the contributions when calculating the skewness 
spectrum. 

In Figure ^ we report both raw and foregrounds- 
cleaned maps data with all corrective terms considered 
compared to C?'^ from Gaussian simulations. 



C. Best Fit Parameter Estimation 

1. Secondary non-Gaussianity only 

At this stage, we make an assumption, relaxed later, 
that there is no primordial non-Gaussianity. 

In Figures || and ^ we plot binned data for each Xi con- 
figuration and for each frequency band with a A/ = 50. 
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FIG. 1: The Icnsing-secondary anisotropy 
cross-correlation power spectrum. We consider two 
secondary effects here with ISW and SZ as these are the 
two dominant lensing-secondary correlations expected. 
The blue solid line is the total contribution when we 
consider both ISW and SZ effects. 



10-1 

10-2 
10-3 
10-4 
10-5 
10-6 
10-7 

10 



XI ISW, 0, V and W band 



100 200 300 400 500 600 




The corresponding error bars are from simulations vari- 
ance and the solid blue line represents a calculation from 
the halo model for the lensing-SZ and lensing-ISW cor- 
relations. These calculations are described in Rcf. |25 . 



In making these estimates, wc also allow for unresolved 
point sources. We calculated the overlap between the 
lensing-SZ and lensing-ISW estimators with the point 
sources given by their shot-noise term of the angular bis- 
pectrum. We parameterize point source amplitude such 
that : 
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' (2/i + l)(2^2 + l) f h h h , 
47r(2/ + l) ^ ' ' 



Bi X 10 sr^ as a function of 



and we consider bgic 
frequency i. 

In order to estimate the relative amplitudes of SZ and 
ISW correlation with the Icnsing potential as well as the 
contaminant contribution from unresolved point sources, 
we compare the data to a model calculation that contains 
lensing bispectra and the overlap between lensing and 
point source bispectra : 
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and Bfi^i^^ is the point source shot-noise bispectrum de- 
scribed in equation (p^. The parameter A denotes the 
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FIG. 2: Corrective terms compared with the raw data 
estimator for Q, V, and W band for Xi with ISW (top) 
and SZ (bottom), respectively. 



total amplitude of the combined lensing-SZ and lensing- 
ISW correlations. 

In doing the above model fit to multi-frequency mea- 
surements of the 2,1 correlator optimized for lensing- 
secondary correlations, we assume that the lensing- 
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FIG. 3: Corrective terms compared with the clean data 
estimator for Q. V, and W band for Xi with ISW (top) 
and with SZ (bottom), respectively. 



FIG. 4: Binned data obtained calculating all the 
corrective terms compared to the skew spectrum of 
simulations for Xi ISW (top) and Xi SZ (bottom). 



secondary cross-correlation amplitude is frequency- 
independent, except for the known frequency dependence 
of the SZ effect taken as part of the model calculation 
of Xi. We assume three different amplitudes for point 
sources in each of Q, V, and W bands of WMAP. 



Our technique as implemented currently does not al- 
low us to separate the lensing-SZ correlation amplitude 
from lensing-ISW amplitude as the two are found to be 
highly degenerate with each other. In future, especially 
with Planck, one may be able to separate lensing-SZ from 
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FIG. 5: Q, V and W band WMAP raw data with error 
bars from simulations variance compared with the 
theoretical models and the best fit results for Xi ISW 
(top) and Xi SZ (bottom). 




FIG. 6: Q, V and W band WMAP clean data with 
error bars from simulations variance compared with the 
theoretical models and the best fit results for Xi ISW 
(top) and Xi SZ (bottom). 
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lensing-ISW based on the SZ frequency dependence. 

To estimate the amphtudes A and Bi we define a 
merit function : 

x'^iy'^ -T-PfC-\y-T-P) (18) 

where y is the data vector (binned Cf'^ for each frequency 
band and Xi configuration), T is the theory matrix and 
P = {A, Bi) is the parameter vector that we want to 
estimate. C is the covariance matrix from simulations. 

To determine the parameters, we minimize the func- 
tion expficitly and obtain the best-fit parameters as : 

P = {T^C-^T)-^{T'^C-^y) (19) 

and the error bars are obtained by the diagonal elements 
of the following matrix : 

a% = [T^C-^ry^ . (20) 

We summarize results related to the amplitude deter- 
mination for different frequency bands in the tables | and 
for raw and foreground-cleaned maps respectively. We 
categorize our results by the choice we make in setting 
Xi while generating C maps. Here, Xi forms a template 
to search for the lensing-secondary cross-correlation. In 
Cf'^, due to weighting, one of the bispectra probes the 
total non-Gaussianity in CMB data that could come from 
a combination of effects from primordial non-Gaussianity 
to lensing and radio point sources. The second bispcc- 
trum forces a certain configuration and the overlap be- 
tween the total and the assumed shape of the bispec- 
trum determines the shape of C^'^ measured from the 
data. By setting a function for Xi, we set the overall 
normalization of the prescribed lensing-secondary cross- 
correlation bispectrum. Thus the amplitude A we de- 
termine from the data is simply the overall amplitude 
of the non-Gaussianity associated with the overlap be- 
tween lensing-secondary cross-correlation and all forms 
of non-Gaussianities that are present in the data. We re- 
move the confusion associated with point sources, which 
is expected to be significant, by explicitly calculating the 
overlap between lensing-(SZ-l-ISW) bispectrum and the 
shot-noise form of the point source bispectrum. 

After accounting for the confusion from point sources 
generated by the overlap of the point source shot-noise 
bispectrum and the lensing-secondary anisotropy cross- 
correlation bispectrum, we find no significant detection 
of the lensing effect in existing WMAP data. We con- 
strain the overall normalization of the lensing-SZ and 
lensing-ISW angular cross-power spectra to be 0.42 ± 
0.86 and 1.19 ± 0.86 in combined V and W-band raw 
and foreground-cleaned maps provided by the WMAP 
team, respectively. The point source amplitude we de- 
termine from the raw map of Q-band with hs-^c = (67.8 ± 
5.4) X 10"^^ sr^ is higher than the estimate by the WMAP 
team with (6.0 ±1.3) x 10"^ /^K^-sr^ H (the value we 
determine is (13.7±1.1) x 10^^ /iK^-sr^ in the same units 
used by the WMAP team). We find similarly a factor of 
2 increase in V-band map as well. 



In the case of clean maps, we find &src = (6.2 ± 
5.4) X 10"^^ sr^, which is smaller than the WMAP 
team's estimate with clean maps for the Q band with 
(4.3 ± 1.3) X 10^5 ^K^-sr^ (the value we determine 
is (1.4 ± 1.1) X 10^^ /iK'^-sr in the same units used by 
the WMAP team). We find similar differences in V and 
W-band as well. 

It is unclear exactly where these differences come from. 
We do not employ the same E-statistic that is optimized 
for point sources as the WMAP team in the present 
study. 

2. Inclusion of primordial non-Gaussianity 

To study the impact of primordial non-Gaussianity we 
now fit the data by modifying equation ( p^ ) to include a 
local form of non-Gaussianity with amplitude /nl : 

D = AxCf^''""~'""-^+B,x£;f'i^^ + /NL>"i'''^"™ , (21) 
where now : 

involves the overlap between lensing-secondary and pri- 
mordial non-Gaussianities with the overall amplitude 
of the primordial non-Gaussianity determined by /nl 

Including primordial non-Gaussianity confuses the de- 
tection of lensing-secondary correlations and leads to a 
factor of 2 degradation in the error of the amplitude of 
lensing-secondary correlation power spectrum (see Ta- 
ble III) . The estimator as constituted is not optimised to 
detect primordial non-Gaussianity, and we find a rather 
weaker limit of /nl = — 13 ± 62 from the clean V-l-W 
maps, with a larger error bar than in the studyof Ref. 
which uses the optimised estimator of Ref. [p4| for pri- 
mordial non-Gaussianity specifically. 

We emphasize that our current study is more focused 
towards a detection of the lensing-secondary correlation 
in WMAP data. In an upcoming paper, we will present a 
combined analysis of three estimators of the 2,1 correla- 
tor optimized separately for primordial non-Gaussianity, 
point sources, and lensing effects. 

In Figure |^ we report the two dimensional countour 
plots showing degeneracies between our best fit parame- 
ters for raw and clean maps from the WMAP team used 
for the data analysis. 

IV. CONCLUSIONS 

We measure the skewness power spectrum of the CMB 
anisotropies optimized for a detection of the secondary 
bispectrum generated by the correlation of the CMB lens- 
ing potential with integrated Sachs- Wolfe effect and the 
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Frequency 


Xi 


A 


Bq 


Bv 


Bw 


xVdof 


Q 


SZ+ISW 


-1.59 ± 1.21 


67.8 + 5.4 






1.67 


V 




0.06 + 1.08 




11.4 + 2.4 




0.85 


w 




1.01 ± 1.06 






5.4 + 2.6 


0.58 


Q+W+V 




0.07 + 0.82 


67.8 + 5.2 


12.4 + 2.2 


5.8 + 2.4 


1.22 


W+V 




0.42 + 0.86 




11.8 + 2.2 


5.2 + 2.6 


0.82 


Q 


ISW 


-1.03 + 1.19 


123.4+12.4 






1.02 


V 




0.33 + 1.06 




20.8 + 6.2 




0.67 


w 




0.99 + 1.05 






10.0 + 7.0 


0.50 


Q+W+V 




0.43 + 0.82 


128.2 + 11.8 


22.6 + 5.8 


11.6 + 6.6 


1.02 


W+V 




0.48 + 0.86 




23.2 + 5.8 


8.4 + 6.8 


0.66 


Q 


sz 


-1.47+1.33 


136.0 + 10.8 






1.68 


V 




0.24 + 1.18 




22.8 + 4.6 




0.84 


w 




1.09 + 1.14 






10.8 + 5.2 


0.60 


Q+W+V 




0.13 + 0.89 


136.2 + 10.4 


25.2 + 4.2 


12.2 + 5.0 


1.23 


W+V 




0.57 + 0.94 




23.8 + 4.4 


10.4 + 5.0 


0.85 



TABLE I: Amplitude parameters estimation using WMAP raw maps. 



Frequency 


Xi 


A 


Bq 


Bv 


Bw 




Q 


SZ+ISW 


2.93 + 1.21 


6.2 + 5.4 






0.55 


V 




0.93 + 1.08 




4.4 + 2.4 




0.70 


W 




1.73 + 1.06 






-1.3 + 2.6 


0.46 


Q+W+V 




1.56 + 0.82 


8.8 + 5.2 


4.2 + 2.2 


-1.4 + 2.4 


0.82 


W+V 




1.19 + 0.86 




5.0 + 2.2 


-1.6 + 2.6 


0.77 


Q 


ISW 


3.32 + 1.19 


17.2 + 12.6 






0.34 


V 




1.16 + 1.06 




5.8 + 6.4 




0.62 


w 




1.81 + 1.05 






-4.3 + 7.0 


0.42 


Q+W+V 




1.76 + 0.82 


25.4+11.8 


5.6 + 5.8 


-5.4 + 6.6 


0.86 


W+V 




1.33 + 0.86 




7.8 + 6.0 


-6.0 + 6.8 


0.67 


Q 


sz 


2.58 + 1.33 


12.2 ± 10.8 






0.62 


V 




0.99 + 1.18 




8.7 + 4.7 




0.70 


w 




1.66 + 1.14 






-2.5 + 5.2 


0.49 


Q+W+V 




1.47 + 0.89 


16.6 + 10.4 


8.2 + 4.4 


-2.2 + 5.0 


0.69 


W+V 




1.22 + 0.93 




10.0 + 4.4 


-3.0 + 5.0 


0.80 



TABLE II: Amplitude parameters estimation using WMAP foreground cleaned maps. 



Sunyaev-Zel'dovich effect. The covariance of our mea- 
surements are generated by Monte-Carlo simulations of 
Gaussian CMB fields with noise properties consistent 
with WMAP. 

When interpreting multi-frequency measurements we 
also take into account the confusion resulting from the 
unresolved radio point sources. We analyze Q, V and W- 
band WMAP 5-year raw and foreground cleaned maps 
using the KQ75 mask out to Imax = 600. 

While with the raw maps we find no evidence for a 
non-zero non-Gaussian signal from the lensing-secondary 
correlation in any of the three bands, we find 2cr and 3cr 



evidence for a non-zero amplitude of both the lensing- 
ISW and lensing-SZ signals in the foreground cleaned 
Q-band map provided by the WMAP team, respectively. 
The point source amplitude at the bispectrum level mea- 
sured with this skewness power spectrum is consistent 
with previous measurements using the optimized skew- 
ness of the WMAP team's analysis and a different form of 
the skewness power spectrum optimized for point sources. 

Finally, as the focus is on secondary non-Gaussianity, 
the estimator is not optimised to detect primordial non- 
Gaussian signals, and we find a limit on local type of 
/nl = —13 ± 62 from cleaned V+W maps. 
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Frequency 


Data 


A 


Bq 


Bv 


Bw 


Inl 


xVdof 


Q 


Raw 


0.39 + 1.99 


69.8 + 5.6 






-95 + 76 


1.68 


V 




-0.51 + 1.78 




11.3 + 2.4 




31 + 70 


0.92 


w 




2.18 + 1.78 






5.6 + 2.6 


-60 + 73 


0.56 


Q+W+V 




1.58 + 1.46 


68.6 + 5.2 


12.2 + 2.1 


5.8 + 2.5 


-70 + 56 


0.60 


W+V 




0.75 + 1.56 




11.9 + 2.2 


5.1 + 2.5 


-16 + 62 


0.85 


Q 


Clean 


2.02 + 1.99 


5.3 + 5.6 






43 + 76 


0.58 


V 




-0.03 + 1.78 




4.2 + 2.4 




48 + 70 


0.73 


w 




3.04 + 1.78 






-1.1 + 2.6 


-67 + 73 


0.42 


Q+W+V 




1.59 + 1.46 


8.9 + 5.2 


4.1 + 2.2 


-1.4 + 2.5 


+ 56 


0.39 


W+V 




1.47 + 1.56 




4.9 + 2.3 


-1.6 + 2.5 


-13 + 62 


0.80 



TABLE III: Amplitude parameters estimation using WMAP raw and clean maps for Xi total and including an extra 

parameter related to primordial non-Gaussianity. 








(a) 



(b) 



FIG. 7: 2-dimcnsional countour plots showing the degeneracies at 68%, 95% and 99.7% confidence levels between 
the best fit parameters from raw (left panel (a)) and clean (right panel (b)) map data analysis for the ISW (top), SZ 

(middle) and joint ISW+SZ (bottom) cases. 
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